Interplay of dendritic avalanches and gradual flux penetration 
in superconducting MgB2 films 



(N 
O 
O 
(N 

o 

Q 



o 
o 

u 

Oh' 



X3 
O 

o 



> 
o 

O 



I 

O 

o 



X 
S3 



D. V. Shantsev^^^ P. E. Goa^, F. L. Barkoyi-^, T. H. Johanseni-*, W. N. Kang, and S. I. Lee^ 
^Department of Physics, University of Oslo, P. O. Box 1048 Blindern, 0316 Oslo, Norway 
^ A. F. loffe Physico- Technical Institute, Polytekhnicheskaya 26, St.Petersburg 194021, Russia 
Institute of Solid State Physics, Chernogolovka, Moscow Region, 142432, Russia 
National Creative Research Initiative Center for Superconductivity, Department of Physics, Pohang University of Science 

and Technology, Pohang 790-784, Republic of Korea 
(February 1, 2008) 

Magneto-optical imaging was used to study a zero-field-cooled MgB2 film at 9.6 K where in a 
slowly increasing field the fiux penetrates by abrupt formation of large dendritic structures. Simul- 
taneously, a gradual flux penetration takes place, eventually covering the dendrites, and a detailed 
analysis of this process is reported. We find an anomalously high gradient of the fiux density across 
a dendrite branch, and a peak value that decreases as the applied field goes up. This unexpected 
behaviour is reproduced by fiux creep simulations based on the non-local field-current relation in 
the perpendicular geometry. The simulations also provide indirect evidence that fiux dendrites 
are formed at an elevated local temperature, consistent with a thermo-magnetic mechanism of the 
instability. 



I. INTRODUCTION 



The dendritic flux instability found in thin films of 
various superconducting materials is a striking but still 
poorly understood phenomenon. It consists in an abrupt 
penetration of magnetic flux into the superconductor 
along narrow branching channels which form irregular 
dendritic patterns on the macroscopic scale. The den- 
dritic instability has been observed by magneto-optical 
(MO) imaging in films of Nb^'^, YBaaCuaOy^-S (induced 
by a laser pulse), and recently in MgB2^~^ and NbaSn^^. 

The instability is believed to be of thermomagnetic ori- 
gin, similarly to the much more explored phenomenon 
of fiux jumping. ^^'^^ Local heating due to flux motion 
reduces the pinning, and facilitates the further motion, 
which may lead to an avalanche process accompanied 
by a substantial temperature rise. This thermal mech- 
anism behind dendrite formation is supported by a re- 
cent experiment, which showed that the instability can 
be suppressed by having a normal metal in contact with 
the superconductor.^ Dendritic patterns of flux® and 
temperature^^ have also been obtained by simulations 
based on the thermal feedback mechanism. 

From the thermo-magnetic nature of this instability, 
one expects that the critical current density jc character- 
izing the flux profile across the dendritic branches reflects 
the elevated temperature at which they were formed. 
Due to a decrease of jc with temperature, these pro- 
files should have a less steep slope than the profiles of 
the regular and smooth penetration from the edges. Sur- 
prisingly, we find that in films of MgB2 the flux profiles 
across the dendritic branches are actually much steeper. 
In this work we investigate this paradox, and study using 
MO imaging the interplay between frozen fiux dendrites 
and the gradually advancing fiux front. 



II. EXPERIMENT 



Films of MgB2 were fabricated on AI2O3 substrates us- 
ing pulsed laser deposition.'^* A 300 nm thick film shaped 
as a square with dimensions 5x5 mm^ was selected for 
the present studies. The sample has a high degree of c- 
axis alignment perpendicular to the plane, and shows a 
sharp superconducting transition at Tc = 39 K. 

The fiux density distribution in the superconducting 
film was visualized using MO imaging based on the Fara- 
day effect in ferrite garnet indicator films. For a recent 
review of the method, see Ref. 15, and a description of our 
setup is found elsewhere. The sample was glued with 
GE varnish to the cold finger of the optical cryostat, and 
a piece of MO indicator covering the sample area was 
placed loosely on top of the MgB2 film. The gray levels 
in the MO images were converted to magnetic field values 
using a calibration curve obtained above Tc- 

The MgB2 film was cooled down to different tempera- 
tures T in zero magnetic field, and then a slowly increas- 
ing perpendicular field, Ba, was applied. For T > 10 K, 
the fiux penetrated the film gradually, and formed the 
critical-state usually found in superconductors with bulk 
pinning. For T < 10 K, the gradual penetration was in- 
terrupted by abrupt invasion of dendritic flux structures. 
Shown in Fig. 1 is a series of MO images visualizing the 
distribution of the perpendicular field B for T — 9.6 K. 
At fields below 14.5 mT, a gradual flux penetration from 
the edges took place. The distribution is then similar 
to the conventional critical-state picture in a perpendic- 
ular geometry:-'^ the fiux concentrates at the edges, seen 
as a bright contour around the sample, and partly pen- 
etrates inwards. The central part and the regions near 
the corners remain flux-free and appear black on the im- 
age. The observed roughness of the flux front is often 
found in MO studies of superconductors and indicate 



1 



the presence of small defects that lead to the fan-like 
flux patterns. At 15 mT a big dendritic flux structure 
abruptly invaded the film (middle panel). The proper- 
ties of such flux dendrites are described in detail in the 
previous studies. ^'''■''"^ Here we emphasize that the exact 
pattern of the dendrite is not reproducible, and thus not 
related to defects. During the subsequent field increase, 
the dendritic structure remains seemingly intact, while 
the flux front continues advancement, like a moving sand 
dune, and by Ba = 36 mT, it covers the dendrite almost 
completely, see Fig. 1 (bottom). 
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FIG. 1. Magneto-optical images of flux distributions in a 
MgB2 film for increasing applied field (only the lower half of 
the film is soon). A dendritic structure appeared abruptly at 
Ba « 15 mT. 

The evolution of the flux density profiles across the film 



during field increase is shown in Fig. 2. The profiles are 
calculated directly from the MO images along the line in- 
dicated in Fig. 1, which is perpendicular to both the film 
edge and two branches of the dendrite tree. The slope 
of the profiles is thus everywhere representing the actual 
|V-B|. The shown profiles cover all the stages of flux pen- 
etration: the gradual advancement of fltix front into the 
virgin fllm (7 and 14.5 mT), the formation of a dendrite 
(15 mT), and the later gradual penetration covering the 
dendrite. One can immediately see that for the dendritic 
branches the profile has an anomalously steep slope. It 
is therefore tempting to conclude that this corresponds 
to an anomalously high critical current density. However, 
an enhanced current density flowing around the branches 
is bewildering. High as it may be during the dendrite 
formation, the current density should relax fast, and not 
exceed jc flowing in the area of regular flux penetration. 
In fact, one would expect that the heating accompanying 
the avalanche reduces the current density. 



Even more surprising is the observed evolution of the 
flux density around a dendrite branch as the flux front is 
approaching and runs it over, see Fig. 3(top). From the 
Figure we again notice the anomalously steep slope of the 
profile across the branch (20 mT curve) compared to the 
slope after the branch has been wiped out (33 and 39 mT 
curves). In addition, one finds that the new fiux com- 
ing to the region from the edge does not simply add up 
to the existing B distribution. Instead, it first destroys 
the existing sharp peak of B{x), so that B is temporar- 
ily decreasing in the vicinity of the branch core. These 
two observations arc in a strong contrast to the behavior 
expected in the usual Bean model, see Fig. 3 (bottom), 
where the slopes of B{x) are fixed, like in a sandpile. 
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FIG. 2. Profiles of flux density for increasing applied field 
obtained from MO images along the line shown in Fig. 1. The 

profile slope across the dendrite branches is much larger than 
that in other regions and, also larger than in the same region 
for Ba = 36 mT. 
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FIG. 3. Top: Profiles of flux density across a dendritic 
branch for increasing applied field Ba-, where x is the distance 
from the film edge. The approaching flux front destroys the 
critical state around the branch leading to temporal decrease 
of local flux density in the peak region. Bottom: Behavior ex- 
pected for a sandpile or the Bean model in parallel geometry, 
which contrasts the experimental results. 



III. SIMULATIONS 



To seek an explanation for the surprising observations 

wc take into account the fact that our sample is a thin 
film. In such a perpendicular geometry the relation be- 
tween B and j is non-local, and this might be responsible 
for the experimental behavior shown in Fig. 3. We ad- 
dress the problem by calculating the evolution of B and 
j distributions for the case of flux front approaching a 
dendritic branch, using flux creep simulations. Focus is 
set on the flux profile development after the dendrite is 
formed, while the dendritic branch itself is introduced 
"by hand" at certain applied field. 

Consider a thin film strip with width 2w along the 
X-axis, and thickness d along the z-axis. We assume 
d <^w and neglect variation of all quantities throughout 
the strip thickness. A magnetic field applied perpendic- 
ular to the strip will then induce a sheet current J = dj, 
and an electric field E directed along the y-axis. From 
the Maxwell equation, one has 



Superconductors in the flux creep regime are usually de- 
scribed by the current- voltage relation. 



E{J,B)=vo\B\ IJ/JepsgnJ, 



(2) 



where n ^ 1, and has the meaning of a vortex depinning 
activation energy divided by kT,^"^ and Vq is the vortex 
velocity at J = Jc- Finally, the nonlocal relation link- 
ing the current and flux density distributions in the strip 
reads 



2- B{x') - Bg jx' - wf 



X — X' 



The simulations start with zero initial conditions, 
B{x,Q) = E{x,0) = J(x,0) = 0, and for t > the 
applied magnetic field is assumed linearly increasing in 
time, Ba = Ba t. To provide correspondence with the 
experiment, the field ramp rate was very slow: Ba <C 
BcVq/w, where Be = I-J-qJc/t^ is a typical values for the 
flux density. The evolution of flux and current density 
distributions J{x,t), B{x,t) is then calculated numeri- 
cally from Eqs. (l)-(3). At the fleld Ba = B* a dendritic 
branch parallel to the film edge is introduced by setting 



B{x) = aB2, xq < X < xo + Ax , 



(4) 



where xq and A give the location and width of the 
branch. Physically, this is equivalent to having the region 
Xq ^ X < Xq + Ax instantly heated to the normal state 
so that it becomes uniformly penetrated by fiux with a 
density proportional to the applied field B*, where the 
coefficient a is determined by the demagnetization fac- 
tor. The simiilations then continue with the same param- 
eters, i.e., we assume that the heated region cools down 
immediately. To mimic our actual experimental situa- 
tion, seen in Fig. 2, the following parameter values were 
chosen: xq/w = 0.5, i.e., the dendrite is formed halfway 
to the film center, at the applied field = 0.58Bc when 
the flux front is located at w 0.75a;o, and to have appro- 
priate height and width of the peak at Xo we set a = 2 
and A/w = 0.03. 

Let us flrst analyze the distributions for the right half 
of the strip, where almost no disturbance is created by the 
dendrite, see Fig. 4. These profiles demonstrate a familiar 
scenario of flux penetration in the perpendicular geome- 
try: nonlinear B{x) with a flux front advancing deeper 
and deeper as Ba increases, and essentially uniform J{x) 
in the flux penetrated area and a considerable current 
also in the Meissner state central part.^^ Such profiles 
are often found in MO studies of thin superconductors,^^ 
and were obtained by flux creep simulations already long 
ago.^° The profiles are also very close to the Bean-model 
result for a thin strip, ^^'^^ which is expected for the large 
n = 25 used in the present simulations. 

Next, we examine the left side of the strip, and how 
the penetration there is perturbed after the appearance 
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of the dendrite. The initial rectangular profile, Eq. (4), at 
B* relaxes very fast because of the high current density 
associated with such an artificial B{x). After this relax- 
ation one finds a triangular shape of the B profile across 
the dendrite - see the thick line in Fig. 4. Note that the 
left slope of the dendrite is steeper than on the right side, 
as found also experimentally, see Fig. 3. As the field con- 
tinues to increase the conventional penetration advances. 



It affects drastically the profile around the dendrite. By 
Ba = 0.75i?c (dashed line), when the flux front has not 
yet reached the dendrite, the peak around .tq is already 
significantly suppressed. As the field reaches 0.92Bc, any 
trace of the dendrite has disappeared, and the flux den- 
sity at xq is only half of its original magnitude. Besides, 
the slope of B{x) around xo became much smaller than 
for the triangular profile. 
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FIG. 4. Flux creep simulations of the evolution of fiux and current density in a thin strip under increasing applied field. At 
the applied field B* = 0.58-Bc a dendritic branch was introduced at xo = 0.5w as described in Eq. (4). The B profiles are in 
qualitative agreement with the experimental results shown in Fig. 2 and Fig. 3. The insets show blow-ups of the area near the 
dendrite. 



IV. DISCUSSION AND SUMMARY 

The present simulations clearly reproduce all main as- 
pects of the experimental behaviour: 

(i) the slope of B{x) across the dendrite is anomalously 
steep; 

(ii) the local flux density in the dendrite core temporarily 
decreases when the flux front approaches; 

(iii) the asymmetry of the triangular profile across the 
dendrite. 

The key to understanding this behaviour lies in the spe- 
cific B — J relation for thin strips expressed by Eq. (3). 
This relation implies an infinite WB where J changes 
abruptly, e.g., near the edge and at the flux front, see 
Fig. 4. A similar situation is present near the dendrite, 
where the current changes direction abruptly. Conse- 
quently, B{x) has anomalously high gradient there too, 
but this is in no way related to having a high local jc- 

Similarly, the non-locality of Eq. (3) is responsible for 
smearing out the B peak at the dendrite. Increasing Ba 
induces additional Meissner currents Jm throughout the 
whole strip, as schematically illustrated in Fig. 5. As a 
result, fiux motion is activated on the side of the dendrite 
where the current density was already jc (the side most 
distant from the flux front). The flux motion proceeds in 
the direction of Lorentz force tending to flatten this side 
of the peak. Earlier MO experiments have revealed sim- 



ilar effects induced by screening currents flowing in the 
Meissner state region of thin films: unexpected flux dy- 
namics around holes^^ and slits^^ ahead of the advancing 
flux front. 




FIG. 5. Schematic of the flux distribution and current flow 
in a thin film where a fiux dendrite is present. 

If the critical state around the dendrite were formed at 
an elevated temperature, the flowing current jc{T) would 
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be correspondingly smaller. As a result, the smearing out 
of the peak will start at a later stage, i.e., when the flux 
front is very close to the dendrite. The magnitude of 
smearing will then also be smaller. In fact, this is what 
we find by comparing the experiments and simulations. 
From the Fig. 4 and Fig. 3 we see that the local decrease 
of B at the dendrite core is « 50% in the simulations, 
but only « 20% in the experiment. This deviation we 
believe results from a local heating in the dendrite area 
during its formation, an effect which was ignored in the 
simulations. 

The work was financially supported by the Norwegian 
Research Council, NorFa and FUNMAT/UiO. 
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